The backbone structure of the Edwards- Anderson spin-glass model 
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We study the ground-state spatial heterogeneities of the Edwards- Anderson spin-glass model with 
both bimodal and Gaussian bond distributions. We characterize these heterogeneities by using a 
general definition of bond rigidity, which allows us to classify the bonds of the system into two sets, 
the backbone and its complement, with very different properties. This generalizes to continuous 
distribution the well known definition of a backbone for discrete bond distributions. By extensive 
numerical simulations we find that the topological structure of the backbone for a given lattice 
dimensionality is very similar for both discrete and continuous bond distributions. Then, we analyze 
how these heterogeneities influence the dynamics at finite temperature and we discuss the possibility 
that a suitable backbone picture can be relevant to describe spin-glass phenomena. 

PACS numbers: 75.10.Nr, 75.40.Gb, 75.40.Mg, 75.50.Lk 



I. INTRODUCTION 

Spin Glasses are the paradigm of systems exhibiting 
both quenched disorder and frustration.^ In such mag- 
netics materials, the static and the dynamical behaviors 
are far from being completely understood. Although in 
the absence of an external magnetic field, experimen- 
tal, theoretical and simulation researches agree on the 
existence of an equilibrium phase transitions at a finite 
temperature, a controversy concerning the true nature 
of the low temperature phase persists. One of theory, 
the mean field picture,^ predicts that spin glasses have 
a non-trivial phase space broken in many ergodic com- 
ponents and with an ultrametric topology. Unlike this 
complex scenario, the phenomenological droplet pictur^ 
postulates a simpler structure for the phase space, with 
only two pure states related to each other by an up-down 
symmetry. 

These two theories have captured the attention of re- 
searchers for at least twenty years, and most of the 
experimental and numerical results have been inter- 
preted in this context. Since the controversy is not 
resolved, other approaches have been proposed to ex- 
plain, within a single framework, many of the results 
reported in the literature4r— More recently, a series of 
studiesii""i^ focusing on the ground-state (GS) topology 
of the Edwards- Anderson ±J spin-glass modeli^ have 
been performed with this same purpose. Based on the 
same spirit as the droplet picture, which focuses on the 
ground state and its excitations, in these works it is as- 
sumed that the spatial heterogeneities of the GS, char- 
acterized by the presence of a backbone, play a funda- 
mental role to describe the low-temperature behavior of 
the system. More specifically, as the GS of the Edwards- 
Anderson ± J model is degenerate, it is possible to define 
the so-called rigid lattice}^' This structure is composed 
by the set of bonds which do not change its condition 
(satisfied or frustrated) in all the configurations of the 
GS. The remaining ones, called flexible bonds, form the 



flexible lattice. Similarly, the sets of spins which main- 
tain their relative orientation in all configurations of the 
GS are known as solidary spins (the remaining spins are 
called non-solidary spins). Both the rigid lattice and the 
solidary spins form the backbone which characterizes the 
GS spatial heterogeneities of the Edwards- Anderson ±J 
model. 

What is observed in these studies is that the back- 
bone structure is closely linked to the static and the dy- 
namical behavior of the Edwards-Anderson ±J model. 
An example of this phenomenon arises in the out-of- 
equilibrium dynamics of the system, which shows strong 
heterogeneities. In particular, the mean flipping time 
probability distribution function has two main peaks cor- 
responding to fast and slow degrees of freedom.— In Refs. 
[ill and |Tj it was shown that both sets are directly re- 
lated to the non-solidary and the solidary spins, respec- 
tively. This suggests that the the backbone structure has 
an influence on the flnite-temperature out-of-equilibrium 
dynamics of the model. In addition, it is important to 
stress that, even below the critical temperature, for long 
simulation times the clusters of non-solidary spins sat- 
isfy the fluctuation-dissipation theorem whereas the sol- 
idary spins violate this relation^iS This suggests that each 
sample could be divided into two structures having very 
different physical properties: the backbone, having the 
properties of a spin glass phase, and the rest of the sam- 
ple in a paramagnetic phase. 

These non-trivial numerical results suggest that a suit- 
able backbone picture can be relevant to describe the 
physics of spin glasses. In Ref. [151 an extensive discussion 
of this idea was carried out. As mentioned there, in order 
to build a more comprehensive theory it is necessary to 
extend the notion of backbone to other disordered and 
frustrated systems. With this purpose, a more general 
deflnition of the rigid structure (see below) of the GS of 
Ising spin-glass models with both degenerate and non- 
degenerate GS was givenji^ As we shall see, a subset of 
this structure could play the role of the backbone of such 



complex systems. 

In order to further investigate these findings, in this 
work we have carried out an extensive numerical study of 
the GS rigid structure of the Edwards- Anderson model, 
for the bimodal ±J and the Gaussian distributions of 
bonds in both the square and the cubic lattices. The 
paper is structured as follows. In Sec. |lT]we present the 
Edwards- Anderson model and the definition of the rigid 
structure. Then, a topological study of this structure 
including a percolation analysis is presented in Sec lIIII 
Finally, the discussion and conclusions are drawn in Sec. 
lYl 



II. THE EDWARDS-ANDERSON SPIN-GLASS 
MODEL AND THE GS RIGID STRUCTURE 

We start by considering the Hamiltonian of the 
paradigmatic Edwards- Anderson spin-glass modelfi^ 
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where the sum runs over the nearest-neighbors sites of ei- 
ther a two-dimensional (2D) square or three-dimensional 
(3D) cubic lattices of linear dimension L and ai — ±1 
are N Ising spin variables. The coupling constants Jij 
are independent random variables chosen from a discrete 
bimodal ± J, 
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distributions, for which the mean value is zero and the 
variance is one. These are the most popular bond distri- 
butions. To avoid confusion, the versions of the Edwards- 
Anderson model where interactions are drawn from Eqs. 
© and (P wiU be called EAB and EAG, respectively. 
Samples (particular realizations of random bond distri- 
bution) in 2D were generated with both periodic-free 
boundary conditions (pfbc) and periodic-periodic bound- 
ary conditions (ppbc), while in 3D only periodic bound- 
ary conditions in all directions were used. 

As discussed in the introduction, for the EAB model it 
is possible to define both the rigid lattice and the clusters 
of solidary spins. These structures can be regarded as the 
backbone of this and other Ising systems with a degen- 
erate GS. Nevertheless, these definitions are not suitable 
to characterize the GS spatial heterogeneities present in 
Ising systems with a non-degenerate fundamental level 
such as the EAG model, because in this case all bonds 
(spins) would be classified as rigid (solidary), forming a 
single homogeneous cluster. 

To understand what changes when a model has a non- 
degenerate GS, following the Ref. ]1^ we begin by con- 
sidering the Edwards- Anderson model with a continuous 



distributions of bonds that consists on the superposition 
of two Gaussian functions of width (variance) e centered 
at J = ±1. We call this the EAB-e model. When e -^ 
the system has the EAB model as a limit case. On the 
other hand, for < e ^ 1 it is reasonable to expect that, 
even though it has a non-degenerate GS (we assume that 
this is the case when bonds are drawn from a continuous 
distribution), the properties of the system, such as the GS 
spatial heterogeneities, will not be very different from the 
corresponding ones for the EAB model. In others words, 
we expect that for a given sample of the EAB-e model 
with a value of e close to but not exactly equal to zero, 
the rigid lattice or the solidary spins of the correspond- 
ing "companion sample" of the EAB model (a sample 
were J — +1 and J — —1 bonds are placed instead of 
ferromagnetic and antiferromagnetic bonds, respectively, 
of the original sample of the EAB-e model) remains the 
backbone of the system. What happens in this particular 
case is that all the GS configurations of the companion 
sample, are distributed in the low-excitation levels of the 
original sample. This suggests that an appropriate back- 
bone for the EAB-e model should be calculated over these 
configurations whose energy is very close to the GS en- 
ergy (in such case the backbone of both models should 
be equal). 

From this particular example we see that, to build a 
general definition of backbone for an arbitrary model, it 
is necessary to consider not only the GS but also the 
low-excitation levels. In particular the "rigidity" of each 
bond should be associated to a parameter taking a con- 
tinuum of values, instead of only two (rigid-flexible) as 
in the EAB model. Such a definition was proposed in 
Ref. [la, and is as follows. Gonsider a sample of the 
Edwards- Anderson model with an arbitrary bond distri- 
bution (discrete or continuous). For each bond Jy we 
define its rigidity as r^ = U*j — U, where U is the GS en- 
ergy of the sample and [/*■ is the lowest energy for which 
the condition of the bond Jij is frustrated (satisfied) if 
it is satisfied (frustrated) in the GS. Then, we will call 
the rigid structure (RS) of a sample, to the lattice where 
each bond Jij has been replaced by their rigidity r^ . As 
shown in the following sections, a subset of the RS com- 
posed by the most rigid bonds seems to behave as the 
backbone of spin-glass models with Ising spins. 
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NUMERICAL RESULTS 



In this section, we present the main topological charac- 
teristics of the RS of the EAB and EAG models in both 
2D and 3D lattices. To calculate this structure, we use 
a similar procedure to that used to determine the rigid 
lattice of the EAB model. ^^'^^ Assuming that one has an 
algorithm for obtaining GS configurations, a scheme of 
the procedure can be described as follows: 

1. For a given sample, a GS configuration C is calcu- 
lated and its energy U is stored. 



TABLE I. Simulations parameters for the 2D EAB and EAG 
models with ppbc (see text). 



TABLE II. Simulations parameters for the 3D EAB and EAG 
models (see text). 
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2. A bond J,-,- is chosen. 



3. The system being in configuration C, one of the ^j^jj „q^^ |-ggg text) 



TABLE III. Parameters for the 2D EAB and EAG models 



spins joined by the bond Jy , i.e. either at or ctj, is 
flipped. This changes the "condition" of the bond 
from satisfied to frustrated or viceversa. 

4. Now the orientation of the spins ai and Cj are 
frozen. 

5. For this "constrained" system a GS configuration 
C* of energy [/*■ is determined. 

6. Then, the rigidity of the selected bond is calculated 

as nj = t/*- - U. 

7. The process is repeated from step 2 until the rigid- 
ity of all bonds have been calculated. 

Notice that the number of GSs required to obtain the RS 
is equal to the number of bonds. 

For this procedure to work one must have an algo- 
rithm that can find any GS configuration of the systems 
involved. But for some systems, if the sample size is not 
very small, only probabilistic algorithms are available, 
i.e. algorithms whose output is a GS configuration with 
a probability smaller than 1. In this case, to improve the 
accuracy of the calculations, the only modification of the 
previous procedure is that, in steps 1 and 5, to obtain a 
reliable GS configuration we perform n independent runs 
of the probabilistic algorithm (evidently, if in step 5 we 
obtain f/* = U, no further runs are performed). 

For lattices with ppbc we use a parallel tempering 
Monte Carlo algorithm j^^i^^ It has recently been shown 
that this technique is a powerful heuristic method for 
reaching the GS of the EAB and the EAG models in 
both 2D and 3D latticesi^ As to determine the RS of 
each sample many independent runs of the parallel tem- 
pering algorithm are needed, we have only been able to 
study lattice sizes of the EAB (EAG) model up to i = 22 
{L = 16) in 2D and L = 10 (L = 8) in 3D. For simplicity, 
in all cases we chose a number of m — 20 replicas of the 
system, and the highest and lowest temperatures are set 
at Ti = 1.6 and Tm = 0.1. In addition, for each model 
with ppbc, Tables HI and HIl list the remaining parameters 
used in our simulations for the different lattice sizes: the 
total number of Monte Garlo sweeps Ngv; , the number of 
sample A'sa and the parameter n. 

On the other hand, for planar lattices it is well known 
that the problem of finding GS configurations can be 
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mapped to a minimum- weighted perfect matching prob- 
lem, which can be solved exactly in polynomial time (i.e 
in time proportional to some power of L).-^ Then, to 
study 2D samples with pfbc, we have used one implemen- 
tation of the Blossom algorithmSS which has allowed us 
to obtain the RS of larger systems sizes. Tables Hill shows 
the corresponding parameters (notice that the Blossom 
is an exact algorithm, so we just take n — 1). In the 
present implementation, that is similar to that used in 
Ref. llSl . many GS configurations must be sampled in 
order to obtain the backbone. Thus, the largest sample 
size that we have studied is L = 60, a smaller size than 
achieved in the cited reference where only the rigid lattice 
was calculated. 

In what follows, first we analyze the size dependence 
of the distribution of rigidity and, by extrapolating, we 
calculate its thermodynamic limit. Next, the percolation 
process displayed for finite lattice sizes is studied in detail 
to infer the structure of the RS of macroscopic samples. 



A. Rigidity distribution 

We begin by analyzing the rigidity distribution, Pi(r), 
for 2D models. Figure [1] (a) shows this distribution for 
the 2D EAB model with ppbc and lattice size L = 22. 
Bonds with rigidity r = are flexible bonds while the 
remaining ones form the rigid lattice. For lattices with 
ppbc, the only possible non-zero rigidity values are r = 4 
and r = 8, i.e., the energy difference between the GS and 
the first and second excited state, respectively. By ex- 
trapolating toward the thermodynamic limit, we obtain 
the asymptotic rigidity distribution, given by: P(0) — 
0.48(1), P(4) = 0.49(1) and F(8) = 0.028(1) [see inset in 
FiglT](a)]. On the other hand, as we can see in Fig. [T](b) 
for samples with pfbc of L = 60, the rigidity distribution 
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FIG. 1. (Color online) Rigidity distributions for the 2D EAB 
model. Samples with (a) ppbc {L = 22) and (b) pfbc {L = 
60). The insets show the peak height of both distributions as 
function of 1/L. 



has two additional peaks at r = 2 and r = 6. These arise 
because excitation energies of AH = 2 and 6 are present 
in lattices with pfbc. In the thermodynamic limit these 
peaks should disappear and both distributions, that of 
systems with ppbc and pfbc, should be equal. Never- 
theless, a linear fitting of the points corresponding to the 
larger sizes (see the inset of FiglT](b)) shows that the dis- 
tribution converges to P(0) = 0.443(1), P(2) = 0.07(3), 
P(4) = 0.451(3), P(6) = 0.009(1) and P(8) = 0.023(1). 
Note that the fractions of bonds with r — 2 and r = 6 are 
non-zero, and the limit values of P(0), P(4) and P(8) are 
close but lower than those obtain for lattices with ppbc. 
These discrepancies are probably due to finite size effects 
caused by the free boundary condition (large finite size 
effects are present in the 2D EAB model with pfbc even 
for larger lattice sizeai^^). We have also found that, on 
average, for each value of the rigidity, approximately half 
of the bonds are ferromagnetic and the other half an- 



FIG. 2. (Color online) Rigidity distributions for the 2D EAG 
model for samples with ppbc (L = 16) and pfbc (L — 60). 
Dots curves correspond to the fits to the Gaussian function 
Eq.(|4]). The inset shows a comparison between Phir) and 
Pl{r*) for the system with pfbc. 



tiferromagnetic. In other words, there is no correlation 
between the rigidity of the bonds and their value. 

We turn now to the 2D EAG model. Figure [5] shows 
a comparison between the rigidity distributions for sys- 
tems with ppbc (L = 16) and with pfbc {L — 60). We 
note that finite size effect are not relevant for this model 
and that Pl (r) is a continuous function taking apprecia- 
ble values within a similar range of rigidity to the cor- 
responding distribution for the 2D EAB model. Again, 
approximately half of the bonds with a rigidity between 
r and r + Ar are ferromagnetic while the other half are 
antiferromagnetic. Thus, we find that for both models 
there is no connection between the sign and the rigidity 
value of the bonds. 

To show that the rigidity distribution is nontrivial (in 
the sense that this distribution cannot be derived from 
any physically reasonable local observable), in the inset 
of Figl5]we show, for the 2D EAG model, a comparison 
between Phii") and the distribution of "local rigidity", 
P^{r*). For each bond Jy of a given sample, we define 
its local rigidity as r*^ = U*^ — U, where as before U is 
the GS energy but now {7* is the smallest of the energies 
of the two configurations obtained by flipping either the 
spin ai or the spin aj (i.e., U*j — U is the minimal value 
of the local excitation needed to change the condition, 
frustrated or satisfied, of bond Jy). Figure [H shows that 
these distributions are very different. In fact, it is easy to 
deduce that for each bond Vij < r* , because the excita- 
tions considered to calculate r, which involve many spins, 
are of lower energy than those considered to determine 



Interestingly the distribution Pl (r) on Figl^] can be 




FIG. 3. Map plot of the "average rigidity lattice" for two 
samples of the 2D (a) EAB and (b) EAG models oi L = 60. 
In both figures the grayscale are the same and the average 
rigidity values are (white), 2, 4, 6 and 8 (black). 



well fitted by a Gaussian function 
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where A and w are two constants. For example, for the 
2D EAB with ppbc we obtain w = 4.89(1) and A = 
1.999(3). 

Another important feature of the RS is its spatial dis- 
tribution. Figure [3] (a) shows a map plot representing 
the RS of a 2D sample of the EAB model with pfbc and 
linear size L = 60. This map plot was generated from an 
"average rigidity lattice" , in which each site value is the 
average rigidity f of the four bonds connecting with this 
site. The grayscale values varying from f = (white) 
to f = 4 (black). We can see that bonds with similar 
rigidity are segregated. Figure [3] (b) shows that this oc- 
curs also for a 2D sample of the EAG model (here we use 
the same grayscale) and that there are small differences 
between their spatial rigidity distributions. 
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FIG. 4. (Color online) (a) Rigidity distribution for the 3D 
EAB model (L — 10). Inset shows the peaks height of the 
distribution as function of 1/L. (b) Rigidity distribution for 
the 3D EAG model for different lattice sizes as indicated, and 
for the 2D EAG model (samples of size L = 16 with ppbc). 
The inset shows a comparison between P{r) and P*(r*) for 
3D samples. 



Next, we analyze the rigidity distribution for the 3D 
models. Figure |4] (a) shows the PL(r) distribution cal- 
culated for samples of the EAB model of size L = 10. 
Now the only possible non-zero rigidity values are r = 4, 
r = 8 and r = 12, which form the rigid lattice (because 
there are few bonds with this latter rigidity value, the 
peak for r = 12 is very small). By extrapolating toward 
the thermodynamic limit we obtain theaymptotic distri- 
bution: P(0) = 0.45(3), P(4) = 0.47(2), P(8) = 0.07(1) 
and P(12) = 0.0016(4) [see inset in Figg](a)]. On the 
other hand, FigH] (b) shows the PL{r) distribution for 
different lattice sizes of the 3D EAG model. Again, we 
observe that the shape of PL{r) is unaffected by finite- 
size effects and that the distribution for L = 8 can be 
well fitted by Eq.ij?]) with parameters w — 5.61(2) and 



A = 1.99(1). A comparison with the rigidity distribution 
of the 2D EAG model, reveals that the range of nonva- 
nishing values of r in 3D is slightly larger, in agreement 
with our observations for the 3D EAB model. Also, by 
comparing Pl(t') with the distribution of local rigidity 
Pl{r*), inset in FigH](b), it is once again clear that both 
functions are very different. Finally, by means of a qual- 
itative analysis of the rigidity structure, we observe that 
a rigidity segregation also emerges in the 3D models. In 
the next subsection, this phenomenon will be evidenced 
indirectly when we study the percolation of the RS of 
such systems. 



B. Percolation of the RS in the EAB model 

In the previous sections we have shown that there are 
small differences between the rigidity distributions of the 
EAB and EAG models, both for 2D and 3D systems. 
Only an increase in the range of r is observed, which is 
a consequence of the change in lattice connectivity. On 
the other hand, we have also determined that the RS is 
not a randomly spatially distributed structure because 
several clusters of bonds with similar rigidity values are 
detected. In this and in the next subsection, we discuss 
how to use the RS to determine which part of the system 
has the relevant properties to be called a backbone. But 
first we need to address the question of which properties 
are to be considered relevant. 

In particular, it seems reasonable to expect that there 
exists a close connection between the topological features 
of the "homogeneous substructures", formed by bonds 
with a similar rigidity, and the equilibrium critical be- 
havior of a given spin-glass model. The main conjecture 
in which this statement is based was proposed in Ref. 
llSl : for a random system with quenched disorder, it is 
the rigidity Vij , and not the corresponding bond strength 
Jij , the quantity that gives the magnitude of the "effec- 
tive interaction" between spins i and j. If this assump- 
tion is true, it is necessary that the backbone of a spin- 
glass model with a finite (zero) critical temperature has 
a percolation cluster with a finite (zero) rigidity value. 
In addition, it is expected that, within this cluster, the 
correlation length diverges at the critical point. In the 
present subsection we study the percolation process dis- 
played within the RS. First, we analyze each one of the 
substructures that formed the RS of the EAB model in 
both, 2D and 3D systems. Then, in the next subsection, 
extending the concept of percolation to lattices with con- 
tinuous links, we perform a similar analysis for the EAG 
model. 

Due mainly to the small systems available to us, which 
have a large dispersion of sizes of the various structures 
we want to analyze, it is very difficult to make a direct 
study of percolation in the EAB model, and to extrap- 
olate it to what could happen in the thermodynamic 
limit,— To overcome this problem, we follow here a more 
complex, but more conclusive, approachJ^i^i^S. In a 
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FIG. 5. (Color online) Percolation probability i?^ for the 2D 
EAB model (ppbc), for the substructures of the RS with (a) 
r = and (b) r = 4. The insets show Phir) and the effective 
thresholds for the percolation criteria / and A, as a function 
of L-i/" for (a) r = and (b) r = 4. 



nutshell, what we do, for each structure to be analyzed, 
is to build a curve of percolation probabilities and to ex- 
tract a percolation threshold from it. Then, we compare 
this threshold with the asymptotic value of the size of the 
structure (corresponding to the thermodynamic limit). If 
this last number is larger than the threshold, we conclude 
that the structure percolates in the thermodynamic limit. 
In the following the procedure is explained in more detail. 
For the set of samples of size L, let us define R^{hr) 
and Rl^{hr)^'^ as the probabilities that a given substruc- 
ture of the RS having a fraction of bonds between hr 
and hr + Ah, percolates along at least one lattice direc- 
tion and percolates simultaneously along all independent 
lattice directions, respectively, and also the arithmetic 
mean of these quantities, Rf;{hr) = [i?^(/ir) + ^i(^r)]/2 
(here the subscript r is the the rigidity of such substruc- 



ture but below it can have a different meaning). For the 
EAB model, we study these percolation probabilities for 
the components in which the RS can be separated, each 
one characterized by a single value of rigidity. In the 2D 
case, for example, the RS of a given 2D sample with ppbc 
may be divided in three clusters of bonds with r = 0, 4 
and 8, whose sizes are Hq, h^ and hs, respectively (these 
quantities are defined as fractions, over the total number 
of bonds). 

Figures[S](a) and (b) show, for the 2D EAB model with 
ppbc, the percolation probability function R"£ of the sub- 
structures with r = and r = 4, respectively, for different 
lattice sizes (for simplicity, we show only the percolation 
criterion A because for this model, it is the quantity less 
sensitive to finite size effects). Calculations were per- 
formed using the algorithm of Hoshen-Kopelman,'^-'^ As 
for small lattice sizes the fractions of rigid bonds such 
as ho and /14 have a very wide distribution)^ the curves 
in these figures extend over almost the entire range of 
this variable, form to 1. In all cases, we have used a 
bin width of Ah — 0.05 and, the corresponding value of 
the percolation probability for samples with a fraction 
of rigid bonds between hr and hr + A/i, was defined as 
the midpoint of the interval, i.e., as hr = hr + Ah/2. 
Error bars were calculated using a bootstrap method^S 
but, from now on, the error bars are omitted for clarity 
when they are equal or smaller than the symbol size. 

The size-dependent behavior displayed for example in 
Fig. \5\ (a), suggests that a bond-percolation process can 
be defined for finite lattice sizes. The crossing point of 
the percolation probability functions allows us to define 
a critical concentration or percolation threshold /iqc- On 
the other hand, for increasing L we have shown in the 
previous sections that the mean fraction of bonds with 
r = tend to a finite value P(0). Combining these two 
results, if hoc < P{0) [hoc > P{0)] we can conclude that 
for a macroscopic sample the substructure analyzed has 
percolated [has not percolated]. 

With this in mind, we analyze the curves in Figs. [S] (a) 
and (b). Although the crossing points are easy to estab- 
lish, to calculate a more precise values for the percolation 
thresholds, we perform for each set a standard analysis 
of the dataJ^i^S First, each curve is fitted with an er- 
ror function using a least-mean-square method. Then, 
the bond concentration at which the slope of the fitting 
curve is largest is taken as an effective threshold /ij^(i), 
where X denotes the percolation criteria used: U, I or 
A. h'^^{L) are expected to follow the law^ 



Kc{L) 



C^'L 



-Ijv 



(5) 



where C^ is a non-universal constant and v is the critical 
exponent associated to the correlation length. 

The inset in Fig. [5] (a) shows, for different sizes, the 
mean fraction of bonds with r = 0, i^L(O), and the ef- 
fective thresholds for the three percolation criteria /, U 
and A as function of L~^l^ . As in a previous workji^ 
here we have used the value oi v = 4/3 of the 2D ran- 
dom percolation;^ In each case, to calculate an estimate 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



FIG. 6. (Color online) The same as Fig. [S] but for the 2D 
EAB model with pfbc. 



of the percolation thresholds we extrapolate towards the 
thermodynamic limit by means of a linear fit. We obtain 
the limits: P(0) = 0.49(1), /i^^ = 0.506(6), h^^ = 0.48(1) 
and ho^ = 0.484(8) (we determine a somewhat differ- 
ent limit of P(0) from the one calculated above, because 
here the fit is carried out using L~'^l^). The closeness 
of these values does not allow us to decide whether the 
substructure with r = percolates, does not percolate 
or is in fact placed exactly at the percolation threshold. 
On the other hand, for the substructure with r = 4, we 
obtain: P(4) = 0.483(6), h\^ = 0.53(1), h\^ = 0.51(2) 
and hf^ = 0.52(1) [see inset in Fig. [S] (b)]. If the per- 
colation threshold h^^ is discarded (in fact this quantity 
is affected by large finite-size effects), we can conclude 
that this substructure has not percolated. Note that the 
substructure with r = 8 has not been considered because 
its size is too small. 

Next we analyze the data for the 2D EAB model with 
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FIG. 7. (Color online) The same as Fig. \S\ but for the 3D 
EAB model. 



pfbc. As in this case the RS of larger lattice sizes are 
available, the situation should be clearer. However, due 
to the existence of large finite size effects caused by the 
free boundary condition, the obtained results lead to a 
similar conclusion. Figures |6] (a) and (b) show the perco- 
lation probability Rl of the substructures with r — and 
r = 4, respectively, for different lattice sizes as indicated. 
In the thermodynamic limit we obtain (see insets in these 
figures): P(0) = 0.45(1), h^^ = 0.49(1), h^^ = 0.46(1) 
and h^^ = 0.48(1); P(4) = 0.46(1), hi^ = 0.59(1), 
h^^ = 0.48(1) and hf^ = 0.53(1). Again, we see that 
the substructure with r = is very close to the percola- 
tion threshold, whereas the other, with r = 4, probably 
has not percolated. 

For the 3D lattices the situation is very different. Fig- 
ures [7| (a) and (b) show for the 3D EAB model, the 
percolation probability i?£ calculated for the substruc- 
tures with r = and r = 4, respectively (for similar rea- 
sons to the above, we omit the study of small clusters 



with rigidity r = 8 and r = 12 because these struc- 
tures clearly do not percolate). In the insets we can 
see a comparison between the corresponding mean frac- 
tions of bonds, and the effective thresholds for only the 
percolation criterion / as a function of L^^''^ (here for 
the 3D random percolation i' — 0.9). In the 3D case, 
for which only small lattice sizes are available, the re- 
maining percolation criteria U and A are not considered 
because these quantities are affected by large finite-size 
effects. In the thermodynamic limit we obtain the fol- 
lowing results: P(0) = 0.45(1) > h^^ = 0.24(1) and 
P(4) = 0.48(2) > hi^ = 0.30(1). As we can see both 
substructures are percolated. Although this cannot hap- 
pen in 2D, for 3D lattices it is possible that two o more 
structures percolate simultaneously. '^4 

It is important to stress that in the EAB models the 
percolation behavior of the substructure of the RS with 
r = 4 is not unlike that found in Ref. [l^ for the rigid 
lattice (recall that the rigid lattice comprises all bond 
with r > 4). This is because bonds with r = 8 in 2D, 
and with r — 8 and r = 12 in 3D, can only form small 
and compact clusters that fill the interstices of the larger 
substructure with r — 4. Then, choosing the backbone 
as being formed by the bonds with r = 4 or r > 4 leads 
to the same percolation and physical properties. 



C. Percolation of the RS in the EAG model 

To perform a similar analysis for the EAG model, first 
it is necessary to extend the above concept of percolation. 
The classical continuum percolation problem has various 
formulations)^ In the 3D problem of voids, for example, 
spherical voids are placed at random in a uniform trans- 
port medium allowing overlapping among them. For val- 
ues of the hole-volume fraction below (above) a critical 
threshold, there (does not) exists an infinite cluster and 
the system is (not) able to support any transport. Here 
the RS of the EAG model constitutes a different sys- 
tem, where the links of a discrete media (the lattice) are 
drawn from a continuous distribution (the rigidity dis- 
tribution). Although many percolation criteria can be 
studied, in this case it is necessary to consider the main 
properties of the RS to define one based on physically 
motivated arguments. 

First, we analyze what happens at finite temperature. 
For a given sample, we define the function Fj(r,T) as 
the mean value of the fraction of bonds with rigidity r 
which, in equilibrium at temperature T, have the same 
condition (satisfied or frustrated) as in the GS. In partic- 
ular, it is easy to see that, when T —>■ 0, those bonds with 
non zero rigidity frozen and then Fj{r, T) — >■ 1, while for 
T —)' oo the mean fraction Fj{r,T) — )• 1/2 (in this case 
each bond is satisfied or frustrated with the same prob- 
ability) . Figure |5] (a) shows the sample average of this 
fraction, F{r, T), for the 2D and the 3D EAB models for 
lattice sizes of P = 16 (with ppbc) and P = 8 respectively 
(bonds with r = does not have a defined condition on 
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FIG. 8. (Color online) F{r, T) as function of r for (a) the 
2D and the 3D EAB models, and for (b) the 2D and the 3D 
EAG models. Curves are given for different temperatures as 
indicated. 



the GS, then we take F(0,r) = 1/2 for any T). Calcu- 
lations for both models were carried out with a parallel 
tempering algorithm, using to = 31 replicas, the highest 
and lowest temperatures set at Tx — 2.0 and Tm = 0.5, 
a number of sweeps of 10^ and the average were taken 
over lO'^ samples of disorder (this set of parameters are 
sufficient to equilibrate both systems). 

With decreasing temperature, we can see in Fig. [8] (a) 
that bonds with high rigidity freeze faster than those 
with low rigidity. This shows that our conjecture about 
the RS is correct: it is the rigidity r^j, and not the bond 
strength Jij , what gives the magnitude of the effective in- 
teraction between spins i and j (it is important to stress 
that in these models all bonds have the same magnitude, 
\Jij\ = 1). In particular for the 3D EAB model and 
close to the critical temperature T^ ~ 1.12;^ we see that 
bonds with r = 12 and 8 are almost completely frozen. 



while only the bonds with r = 4 are affected by thermal 
fluctuations. But, as was shown above, it is precisely 
this latter substructure what percolates in the thermo- 
dynamic limit. Then, it is reasonable to assume that, at 
the critical point, the correlation length diverges mainly 
throughout the region formed by the bonds with rigid- 
ity r = 4 (marginally, bonds with a greater rigidity than 
4 could connect the percolation cluster to the islands of 
finite size, all belonging to the same substructure with 
r = 4). 

On the other hand, in Fig. [8] (b) we show the same 
function F{r,T) for the 2D and the 3D EAG models, 
for lattice sizes of L — 16 and L = 8 respectively. To 
make the calculations, we have used a parallel tempering 
algorithm that was run with the same parameters as for 
the EAB models. We observe a similar behavior to that 
shown in Fig. [S] (a), but now for models with continuous 
bond and rigidity distributions. Notice that close to Tc ~ 
0.95,'^^ bonds with a rigidity between r « 2 and r sa 4 
are affected by thermal fluctuations. 

So far, we can make two observations that are relevant 
to our present study. On the one hand, as mentioned at 
the end of the previous subsection, although the proper- 
ties of each subset of the RS are important, for practical 
purposes it is enough to define the backbone as being 
formed by the bonds with r > 4. On the other hand, as 
shown in FigslS] (a) and (b), at the critical temperature 
the regions of high rigidity are almost completely frozen. 
Therefore, whether bonds of high rigidity are included or 
not in the backbone, we expect that both choices lead to 
the same percolation and physical properties. 

The arguments given in the previous paragraph sug- 
gest that it would suffice to break the RS into two 
sets, one formed by the bonds with rigidity r > r^^^, 
^(''min)' ^-'^'^ another set ^*(''niin) comprising the re- 
maining bonds, where r-^^-^ is a given rigidity value. Note 
that this procedure is similar to separate the system into 
two lattices, fiexible and rigid. Although for the EAB 
models this approach does not provide new information, 
for the EAG models it supplies a better way to treat the 
percolation problem of the RS (given the small lattice 
sizes, it is not possible to analyze the percolation prop- 
erties in the EAG model by dividing the RS into many 
substructures). 

Then, to study the percolation of the RS of the EAG 
models we proceed as follows. For each sample of linear 
size L and for a given value of r^^, the substructure 
^(''min) contains a particular fraction of bonds h^. This 
quantity has a wide disorder distribution but its average 
converges to a value given by 



^(''min)~^i(^min) 



PUr) dr, 



(6) 



which does not depend appreciably on the size. Similarly 
to what we have done for the EAB models, for each of 
the possible values of T'jnin' ^^ calculate a percolation 
threshold ^a;c(''niin) ^'^'^ then we compare this number 
to the corresponding -'^('"niin)' ^^ ^^^ ^ -^ \^^<^ ^ -^^ 
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FIG. 9. (Color online) Percolation probability i?^ curves for 
(a) the 2D (ppbc) and (b) the 3D EAG models, for r^^^ = 2.5 
and fjjjjjj = 4 respectively. The insets show for each model, 
the comparison between hxc and X, and between /i*^ and 
1 ~ X, as functions of r™;„. 



the set ri [does not] percolate. Figure [5] (a) shows as 
an example, the curves of the percolation probability R^ 
for the 2D EAG model with ppbc, for r^^-^ = 2.5. In 
this case the percolation threshold is /ia;c(2.5) = 0.48(2) 
(calculated as before, through a standard analysis of the 
data) and X(2.5) = 0.31(1), therefore the set f}(2.5) does 
nor percolate. For both low and high values of the rj^ij^, 
it is difhcult to calculate a reliable percolation threshold, 
because the set f^(''min) for most finite-size samples al- 
ways percolate or does not percolate, respectively. On 
the other hand, for the lattices with pfbc, it also is com- 
plicated to obtain this quantity because the curves of the 
percolation probability (for any criterion) are very steep 
{hx has a narrow distribution for the biggest lattice sizes). 
The inset in Fig. |9] (a) shows the comparison between 
hxc and X as function of r^^. Extrapolating we can 



sistent with a Tc = for this model (this picture is similar 
to that found for the 2D EAB model) . Note that h^c in- 
creases when T-jj^jjj decreases. In order to understand this 
behavior, let us consider two different values of '"niin' ''i 
and r2, with ri > r2 but close to each other. If for a given 
sample the set r2(ri) percolates and the corresponding 
hx ^ hxc{ri) then, although r^^^^^ is decreased to r2, this 
set will still percolate. Also, for another sample where the 
set ri(ri) does not percolate and hx <C hxci^i), we ex- 
pect that rj(r2) does not percolate. On the other hand, 
for samples which have a value of hx close to hxcifi), our 
data show that in most cases the percolation state (per- 
colating or not percolating) of both sets ri(ri) y Sl(r2) 
does not change. This is probably because in these situ- 
ations an increase in the number of bonds composing the 
set n does not appreciably alter the size of the largest 
cluster (which we expect to be fractal), but serves to in- 
crease the size or the number of the smaller islands. In 
all cases, when we reduce the value of r^j^j^ the frac- 
tion hx is increased [the set Q{r2) has more bonds than 
the set ^(ri)] and the curves of the percolation proba- 
bility move to the right (to higher values of hx) without 
changing appreciably its shape. As a consequence, the 
percolation threshold increases: hxc{ri) < hxc(T2)- 

Now, let us consider the complementary set ^*(''niin) 
composed by the bonds with r < r^^. For this set it is 
possible to calculate a percolation threshold ^i^c('"min)- 
Comparing this quantity with the fraction 1 — ^('"niin)' 
we can find whether the bonds with the lowest rigidity 
percolate in the thermodynamic limit. The inset in Fig.[S] 
(a) shows this comparison. As we can see for the 2D EAG 
model, when r-j^jj^ > 1.3 the set ^*(^niin) percolates. 

For the 3D EAG model we obtain a different behav- 
ior. In Fig. E] (b) we can observe for the set ^(j^uiin)' 
that the percolation probability i?£^ curves for r-^^ — 4 
intersect at hxc{4:) = 0.16(1), a value very close to 
X{4) = 0.163(1). In turn, the inset in this figure shows 
that hxc and X cross at a value of r^^ close to 4. On 
the other hand, for the complementary set ^*(?'niin)' ^^^ 
inset in Fig. [5] (b) also shows that the curves of /i*^ and 



1 ~ X cross at r,. 



0.6. 



D. The backbone of the EAG models 

In the previous sections we have shown that the per- 
colation properties of the RS for the EAG and the EAB 
models are very similar. In 2D, seemingly there is only a 
percolation set that contains the bonds with the lowest 
values of rigidity. On the other hand, for 3D systems 
there is either a particular value or a range of r^^ for 
which both sets f^C^^min) and f^*(''niin) percolate. This 
suggests that by using a suitable value (or range of val- 
ues) for Tjjjjjj the sets ^(j^niin) ^^'^ ^*(^min) '^^^ have 
similar properties to those of the rigid and the flexible 
lattices, respectively, of the EAB models. In this way, 



11 



■b~u \ 


0.30- 


Q ; Backbone 




■ \ 


0.25- 


0.20- 
0.15- 


\ : V 2DEAG(ppbc) 

\ \ \ ^-J^w 


0.10- 


0.05- 
0.00- 


A ^O 

A^A-A.A A A A A A R-Cj3-Kl-m 



4 
r 



0.35- 

0.30- 

0.25- 

0.20- 

0.15 -p^ 

0.10 

0.05- 

0.00 



Backbone 



.o-°-°--o \ 



A 



A : 






3DEAG 

-n-P,{r) 

-^-P[{r) 



^A 



"A 



"A, 



'"A, 



~A- 



(b). 



-D-, 



□ -r 



FIG. 10. (Color online) Rigidity distributions PL(r), Pfir) 
and P[(r) (see text) for (a) the 2D (L — 16 with ppbc) and 
(b) the 3D (L = 8) EAG models. The vertical dotted lines 



marks the value of r-^^j^ = 1.3. 



we could define the backbone of these models as the set 
^(''min)- ^^ ^^^ following we show how a suitable value 
(or range of values) for r^^^jj-^ can be found. 

In Figs. [10] (a) and (b) we can see, respectively, for the 
2D and the 3D EAG models, a comparison between the 
rigidity distributions for all bonds, Pl(?'), for the bonds 
that are satisfied in the GS, Pf{r), and for the remaining 
bonds that are frustrated in the GS, Pf{r). Obviously 
the distributions must satisfy PLir) — -?/(»') + Pl {''')■ 
Note that the distributions P£ (r) and P[ (r) are not nor- 
malized to unity because their integrals are equal to, re- 
spectively, the fractions of satisfied and frustrated bonds 
of the GS. In both cases, for high values of r^^^, al- 
most all the bonds in the set ^(''jxiin) ^-'^^ satisfied in 
the GS. However, with decreasing r^^^^^^, although more 
and more satisfied and frustrated bonds are incorporated 
into the set ^^(''jxiin)' clearly this structure remains with 



a low frustration. A change of trend occurs when the 
Pf{r) distribution reaches its maximum: for lower val- 
ues of ^j^in ^^"^ ^^^ ^(^min) begins to incorporate highly 
frustrated regions. Remarkably, this maximum happens 
approximately at r^^ « 1.3 for both, the 2D and the 
3D EAG models. 

As we have seen above, for the 2D EAG model this 
behavior is significant since at r^^^y^ = 1.3 the comple- 
mentary set il*(1.3) is very close to (or is at) the per- 
colation thresholds, while the set $7(1.3) does not perco- 
late. This picture occurs at a precise value of rigidity and 
corresponds to that observed in the 2D EAB model for 
''min = 4- Therefore, as suggested in the Fig. [10] (a), it 
is possible to define the backbone of the 2D EAG model 
simply as the set £7(1.3). This choice produces a back- 
bone having the 60% of the bonds of the system and a 
fraction of frustrated bonds in the GS of only 0.085. 

On the other hand, for the 3D EAG model the prob- 
lem of determining a single value of ^'u^jj^ is not so sim- 
ple. Nevertheless, as we show below, choosing again 
''min ~ 1.3 we obtain a structure with the appropriate 
characteristics to become a suitable backbone. Note that 
in this case this could also be achieved by choosing r^^^ 
within the range of values that goes from 0.6 to 4. 

In the EAB model, several observables behave differ- 
ently when they are calculated within or outside of the 
backbone. Such is the case of the average energy per 
bond. In Ref. [l5| it was shown that, for example, when 
this quantity is evaluated outside the backbone (in the 
flexible lattice) it has a minimum at low temperature (see 
Figs. 17 and 18 in this work). To perform a similar anal- 
ysis for the EAG models, let us define u(T), uq{T) and 
Ufi" (T) as the average energies per bond at temperature 
T of, respectively, the whole system and the sets il and 
n* for a given r^jy^y^. 

Figure [TT] (a) shows, for the 2D EAG model, these 
energies as function of T for r^^^^^ — 1.3. Just as in the 2D 
EAB model,J^ we observe a very broad minimum in the 
curve of Ufi-iT). For higher (lower) values of r^^^^, the 
inset in that figure shows that the minimum disappears 
(becomes narrower). Again, we confirm that choosing 
^min — 1-3 produces a suitable backbone for this model. 

However, for the 3D EAG model, Fig. [11] (b) and its 
inset show that for some values of r^^yy^ around 1.3, the 
curves of uq* (T) display a minimum. Our calculations 
suggest that a suitable backbone could be defined in the 
range [0.6 — 2.0], since all these structures have similar 
topological characteristics. However, to simplify subse- 
quent studies, we choose r^j^yj^ = 1.3 to define the back- 
bone of the 3D EAG model. The backbone thus defined 
has 64% of the bonds of the system and a fraction of 
frustrated bonds in the GS of 0.163. 

Note that the percentage of bonds comprising the 
backbone is only slightly larger in 3D than in 2D, while 
in 3D the fraction of those bonds that are frustrated is 
twice that in 2D. The same is obtained for the EAB mod- 
els: the backbone (henceforth, for the EAB models we 
will assume that the backbone consists of all bonds with 
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FIG. 11. (Color online) Average energies per bond u{T), 
uniT) and un*{T) (see text) for (a) the 2D (L = 16 with 
ppbc) and (b) the 3D {L — 8) EAG models. In both cases 
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rigidity r > 4) comprises approximately 52% of the sys- 
tem in 2D and 54% in 3D, and the fraction of frustrated 
bonds is 0.05 in 2D and 0.1 in 3D.l^ 

Finally, we analyze the internal structure of the back- 
bone studying the cluster number distribution, ng., i.e. 
the number of clusters of size s. For the random bond 
percolation at the critical concentration, it is expected 
that this distribution follows a power law 



Us oc s 



(7) 



where r is a critical exponent i^ Because large samples 
are needed, in 2D we have restricted our analysis to lat- 
tices with pfbc of size L = 60. Figure [T^ (a) shows, for 
the 2D EAB model, the cluster number distribution cal- 
culated for a range of hx centered at 0.55, the mean frac- 
tion of bonds with rigidity r > 4. Fitting the curve we 



FIG. 12. (Color online) Cluster number distributions for the 
EAB and EAG models for (a) 2D (L = 60 with pfbc) and 
(b) 3D (L = 8) lattices. The distributions are calculated for 
different ranges as indicated. 



obtain t = 1.95(5). On the other hand, we have also cal- 
culated this distribution for the backbone of the 2D EAG 
model (rj^in = 1-3), for two different ranges of hx- one 
centered at /i*j,(1.3) = 0.68 and another at X(1.3) = 0.6 
[see inset in Fig. ^ (a)]. For these ranges we obtain, re- 
spectively, T = 2.03(8) and r = 1.98(4). Whatever the 
model or the range, the exponent values are very close to 
T — 187/91 « 2.05, the corresponding exponent for the 
2D random percolation universality class.'^'^ 

For 3D systems Fig. [T^] (b) shows the cluster number 
distributions for both the 3D EAB model for a range cen- 
tered at 0.54 (the mean fraction of bonds with rigidity 
r > 4), and the 3D EAG model for a range centered at 
X(1.3) = 0.64 [see inset in Fig. M (b)]. We obtain, re- 
spectively, T — 1.71(9) and r — 1.71(5), values that are 
clearly different from the accepted exponent r — 2.2 of 
random percolation in 3D;^ This difference is presum- 
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ably due to the fact that, unhke the 2D case, here the 
mean size of the backbone tends to a value far from the 
percolation threshold. Nevertheless, for the small sizes 
considered, we see that the cluster number distributions 
for both 3D models, follow a power law with the same ex- 
ponent. These results show that the internal structure of 
the backbones of the systems with bimodal and Gaussian 
bond distributions are similar. 



IV. DISCUSSION AND CONCLUSIONS 

In the present work we have studied the topology of 
the rigid structure or RS, which is defined in the GS of 
the Edwards- Anderson spin-glass model or, in general, in 
the GS of all Ising-like systems. So far, the main results 
are that this structure allows us to find a suitable back- 
bone for both discrete and continuous bond distributions 
and, for a given spatial dimension, that these backbones 
share the same topological properties. In the following 
we discuss in more details two additional and important 
topics. 

In the previous sections we have shown that the con- 
nection between the bond strength J^ and the rigidity 
rij is not trivial. In principle, this means that given a 
particular bond it is not possible to predict its rigidity 
value. As we have seen, for the bimodal distribution this 
is true. To further analyze this issue for the Gaussian dis- 
tribution, we calculate for the 2D and the 3D EAG mod- 
els the probability distribution that a bond has strength 
J and rigidity r. Dividing this function by Dq(J), the 
Gaussian bond distribution, we obtain the conditional 
probability density W{J,r) (the conditional probability 
is chosen because the bond distribution is not uniform 
and then bonds with a low strength are more numerous 
than those with high strength). 

Figures [13] (a) and (b) show the map plots of W{J,r) 
for, respectively, the 2D and the 3D EAG models. First, 
one can see that for high values of |J| and r both dis- 
tributions develop two arms. What this shows is that 
for most of these bonds their strength is proportional to 
their rigidity. This is to be expected, since these bonds 
must be surrounded by others of lower strength (which 
are more numerous) and then we expect that the changes 
in the GS energy (the rigidity) , produced by changing the 
condition of a bond of great magnitude, should depend 
primarily on the value of | J|. In other words, the rigidity 
of a high-strength bond seems to be, on average, a trivial 
quantity)^ 

On the other hand, for intermediate and low values of 
\J\ and r, the distribution of W{J,r) for both the 2D 
EAG and the 3D EAG models, has a "square structure" 
in which the proportionality between the strength of a 
given bond and its rigidity is lost. For example in the 
3D case. Fig. [T3] (b) shows that W{J, r) is almost con- 
stant on the region where bonds lie in the approximate 
ranges — 1 < J < 1 and 1 ;$ r < 2. The same applies to 
bonds with | J| sa 1 which have, with equal probability, a 




FIG. 13. Map plot of the conditional probability density 
W{J,r), for (a) the 2D (L = 60 with pfbc) and (b) the 3D 
{L — 8) EAG models. The horizontal dotted lines marks the 



value of fjjjjjj = 1.3. 



rigidity between < r < 2. These examples show that 
the intermediate region of parameters, where many of the 
bonds that make up the backbone are located, and where 
the thermal critical (in 3D) fiuctuations are important, 
is non-trivial. 

The other topic is concerned directly with the back- 
bone picture. As discussed in the introduction, this phe- 
nomenological theory has been mainly inspired in previ- 
ous numerical studies of the EAB models, in which there 
has been shown that there is a connection between the 
GS spatial heterogeneities (characterized by the back- 
bone) and the finite-temperature behavior of these mod- 
els. The main conclusion of these works is that the back- 
bone plays a fundamental role in the physics of some dis- 
ordered and frustrated systems. Although the RS stud- 
ied here and the corresponding backbone obtained from 
it still remains as structures that are defined in an ar- 
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bitrary way, some steps can be taken in order to clarify 
their physical interpretation. 

Previous papers have shown that the backbone and 
its complement seem to influence both the equilib- 
rium and the out-of-equilibrium dynamics of the EAB 
modelsJ^ii^"— Here, Figs. [11] (a) and (b) show that the 
same happens in the EAG models (in the equilibrium dy- 
namics). The problem with this interpretation is that it 
might be considered odd that the backbone, or the RS 
obtained by analyzing only the GS and the first-excited 
energy levels, may influence the dynamics which takes 
place at finite temperature, exploring configurations cor- 
responding to highly excited energy levels. To address 
this point, another interpretation of the RS was proposed 
above: for a random system with quenched disorder, the 
rigidity r^j and not the corresponding bond strength Jij , 
is the quantity that give the magnitude of the effective 
interaction between spins Ui and Uj. 

This latter point of view is motivated by the following 
arguments. The backbone picture should be valid for a 
broad class of the systems, ordered as well as disordered. 
A typical ordered system is the 2D ferromagnetic Ising 
model, whose Hamiltonian is given by Eq. ([1} but with 
Jij — J > 0. Its RS is trivial to obtain: it is easy to 
see that all bonds have a rigidity of r = 8 J which is 
proportional to J. Those phenomena occurring at finite 
temperature, as the continuous phase transition whose 
critical temperature is Tc ~ 2.269 J, also depend on the 
bond strength. 

We clearly see that the RS of the ferromagnetic Ising 
model does not project at finite temperature and, as the 
rigidity is proportional to the bond strength, both quan- 
tities are a measure of the interaction between spins ai 
and aj. Simply the protocol used to calculate the RS, 
enables us to determine all effective couplings. Logically, 



the same should happen in a spin glass where we expect 
that the quenched disorder and the frustration eventu- 
ally weaken both, the bond rigidities and the effective 
interactions between all pairs of spins. 

Finally, we mentioned that the topological character- 
istics of the backbone are also relevant to understand- 
ing tjie physical behavior of other random systems. In 
Ref. i37i, for example, Tsomokos et al. predict that, if the 
backbone of the 2D EAB model does not percolate, in the 
random-field toric code model there exists a new inter- 
mediate quantum phase where topological and spin-glass 
order coexist. Our calculations suggest that this could 
also be the case for both 2D models. In addition, in a re- 
cent study of the out-of-equilibrium dynamics of the 2D 
± J Potts model at low-temperature, numerical evidence 
has been found that hints at the existence of an underly- 
ing backbone structure for this systemi^ Unfortunately, 
the RS studied here is defined only for models with Ising 
spins (the q = 2 case in the ±J Potts model correspond 
to the EAB model). Although a general procedure to 
obtain the backbone of an arbitrary system has not yet 
been found, we believe that the progress made here rep- 
resents a significant step in this direction. In addition, we 
speculate that a further generalization of the rigid struc- 
ture analyzed here will allow us to identify the backbone 
of more complex disordered models. 
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